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We consider a parametrization of the lepton mixing matrix in which the 

CN ■ deviations from maximal atmospheric mixing and vanishing reactor mixing are 

\f^ • obtained in terms of small corrections from the charged lepton sector. Rela- 

^^ ' tively large deviations for the reactor mixing angle from zero as indicated by 

T2K experiment can be obtained in this parametrization. We are able to fur- 
r^ I ther reduce the number of complex phases, thus, simplifying the analysis. In 

addition, we have obtained the sides of unitarity triangles and the vacuum os- 
cillation probabilities in this parametrization. The Jarlskog rephasing invariant 
. , . measure of CP violation at the leading order has a single phase difference which 

\^ . can be identified as Dirac-type CP violating phase in this parametrization. 

1 Introduction 

Results from a variety of solar, atmospheric and terrestrial neutrino oscillation exper- 
iments [1] have constrained the form of the lepton mixing matrix U [2]. The lepton 
mixing matrix is given by 

U = UJU. (1) 
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where Ui and 11^ are both 3x3 unitary matrices such that Ui arises from the diago- 
nahzation of the charged lepton mass matrix (Mi) while 11^ diagonahzes the neutrino 
mass matrix (M,,). For three lepton generations, the 3x3 unitary matrix U in Particle 
Data Group (PDG) [3] parametrization is given by 
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where S13 = sin ^13, C13 = cos ^13 with ^13 being the reactor angle, S12 = sin ^12, C12 
= cos ^12 with 9i2 being the solar angle, S23 = sin ^23, C23 = cos ^23 with 6'23 being 
the atmospheric mixing angle and 6 is the Dirac-type CP violating phase. The phase 
matrix p = diag(l, e*"^/^, e*"^/^) contains the Majorana-type CP violating phases ai 
and 0:2 which do not affect neutrino oscillations and are not directly accessible to 
experimental scrutiny at present. Current data is consistent with the tribimaximal 
(TBM) mixing [4] 
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which has been derived using family symmetries [5]. In addition to TBM, there 
are other mixing schemes which can reproduce the observed leptonic mixing pattern 
which include the two Golden Ratio (GR) mixing schemes where the mixing angles 
are for GRl: ^12 = tan'^(l/v9), ^23 = 7r/4, ^13 = [6], GR2: ^12 = cos-i(v9/2), ^23 
= 7r/4, ^13 = [7] where y? = (1 + 75)72, Hexagonal Mixing (HM): 6^2 = n/Q, 623 
= 7r/4, ^13 = [8], Bimaximal Mixing (BM): ^12 = 7r/4, ^23 = vr/4, ^13 = [9]. All 
these mixing schemes can arise from mass independent textures also known as form 
diagonalizable textures [10] and lead to a predictive neutrino mass matrix structure 
which contains just five parameters (the three neutrino masses and two Majorana 
phases). All the above mixing scenarios have the same predictions for the reactor and 
atmospheric mixing angles viz 6'i3 = and ^23 = 7r/4, whereas their predictions for 
the solar mixing angle ^12 are different. Thus, the above mixing matrices are common 
upto a mixing matrix 
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arising from a mu-tau symmetric neutrino mass matrix. It is highly unlikely that 
any of the above mixing schemes is exact since there are already hints of a non-zero 



reactor mixing angle ^13 [11]. Recently, the T2K collaboration has observed possible 
indications of the u^ — ?> z/g appearance and reported the following ranges for ^13 [12] 

5.0° < ^13 < lQ.O°forNH (5) 

5.8" < ^13 < 17.8" forlH (6) 

at 90% C.L.. Moreover, the best fit value of ^13 is found to be ^13 ^ 9.7° for NH and 
^13 ^ 11° of IH, thus, implying large deviations from ^13 = 0° in the above mentioned 
mixing scenarios. Therefore, it becomes important to develop a parametrization of the 
lepton mixing matrix in which such deviations are manifest. A natural possibility to 
obtain a phenomenologically viable neutrino mixing matrix and to generate non-zero 
^13 and non-maximal 623 is to assume that these deviations come from the charged 
lepton sector. Such an assumption has been made earlier to generate deviations from 
bimaximal mixing [13, 14] and tribimaximal mixing [15, 16, 17, 18]. 

2 Formalism 

A general 3x3 matrix contains 3 moduli and 6 phases [19] and can be represented as 

U = e'^PUQ (7) 

where P = diag(l, e**^^, 6**^^) and Q = diag(l, e*''^, e*''^) are diagonal phase matrices 
having two phases each, U is the matrix containing 3 angles and one phase and 
has the form of U (except for the phase matrix p) in Eq.(2). In general, when 
charged leptons also contribute to the mixing, the lepton mixing matrix contains 6 
real parameters and six phases [14]. As pointed out earlier the two Majorana phases 
are unlikely to be measured in the present and the forthcoming experiments, so that 
we may dispense with the Majorana phases at least for the present by considering 
the Hermitian products MiMJ and M^Ml. Here, two points are in order: 
1.) Since we are considering mass independent textures of the neutrino mass matrix, 
thus, Mj, and M^Ml are diagonalized by the same diagonalizing matrix, so we can 
consider the product M^M^. 

2.) The deviations of the charged lepton mass matrix from diagonal matrix are in 
any case considered to be arbitrary so the choice of product MiMJ can be made. 
By using these Hermitian products we not only dispense with the unnecessary burden 
of Majorana phases but we are also able to remove one additional phase from the 



lepton mixing matrix, thus, simplifying the subsequent analysis. The lepton mass 
matrices can be diagonalized as 

Ml = UiMfUJ, (8) 

M, = U.MtU^. (9) 

Thus, the product MiMJ becomes 

MiMJ = UiMfuiURMfu} (10) 

= UiiMffU} 

which can be written as 

MiM} = e"^'PiUiQi{MtfQ\u!pU-"^' (11) 

= PiUi{MffUi^ P} 

using Eq.(7). Similarly, for the product M^^Ml we obtain 

M^Ml = P,U,{MtfuJ Pi (12) 

We can absorb two phases from Pi and one phase from P^ in the left handed lepton 
fields and the resulting lepton mixing matrix is given by 

U = UiPuUu (13) 

where Ui and U^ contain three real parameters and one phase each while P^ contains 
one phase Py = diag(l, 1, e*"^). Thus, in this formalism U is expressed in terms of six 
real parameters and three phases. 

In the present work, we discuss a parametrization of the lepton mixing matrix which 
allows for the large deviations from ^13 = and has the form of Eq.(4) at zeroth 
order. Here 9[2 can have the values sin^^(l/v^) for TBM mixing, tan~^(l/(/9) for GRl 
mixing, cos^^{ip/2) for GR2 mixing where y? = (1 + v5)/2, 7r/6 for hexagonal mixing 
and 7r/4 for bimaximal mixing. Deviations from the above mentioned scenarios are 
parametrized in terms of charged lepton corrections represented by small parameters 
having magnitude of the order of Wolfenstein parameter A ~ 0.227 or less. In the 
small angle approximation, we have to the first order in eij 

sineij ~ eij, cose^j ~ 1 eij < 0.227 (14) 

where i,j = 1,2,3 and i < j. 



3 Deviations from exact mixing schemes 



In this section, we obtain expressions for lepton mixing observables in terms of the 
charged lepton corrections. The neutrino mixing matrix Ui, is assumed to have the 
form of Eq. (4) except for the phase matrix p. The charged lepton mixing matrix to 
first order ^ in terms of small parameters is given by 
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and the resulting lepton mixing matrix has the form 



U = Ui^PM, 
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The lepton mixing angles are related to the elements of the mixing matrix as 
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The mixing angles in this parametrization, to first order in small parameters are given 

by 
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It can be seen that the reactor mixing angle along with the atmospheric mixing an- 
gle is independent of 6[2 and the deviation of the atmospheric mixing angle from 
maximality depends only on the small parameter €23 in the first order corrections. 
We restrict the ranges of perturbation parameters by using the recent global analysis 



^second order corrections to the mixing matrix elements and other relevant quantities have been 
discussed in the Appendix 



[20] which incorporates the T2K [12] and MINOS [21] results. Allowed numerical 
ranges for the perturbation parameters ei3 and £23 at 3cr for all mixing scenarios are 
—0.22 < (ei3, 623) < 0.22. The range of eu at Scr is 

-0.20 < ei2 < 0.17 TBM 
-0.165 < ei2 < 0.22 GRl 

-0.22 < ei2 < 0.17 GR2 

-0.15 < ei2 < 0.22 HM 

-0.22 < ei2 < BM. (19) 

The Jarskog CP violation rephasing invariant [22] is given by 

sin 2^;2ei3 sin(5i3 - 0) 
JcP = ^ . (20) 

An important point to note here is that the Jarskog rephasing invariant Jcp to the 
first order contains a single phase difference which is the relevant Dirac-type CP 
violating phase in this case. However, it does not necessarily coincide with the Dirac- 
type phase in the standard parametrization as pointed out in Ref. [17]. Also, in Ref. 
[16, 17], the expression for Jcp contains two different phases in the leading order 
term and further assuming CKM like hierarchy among the perturbation parameters, 
the relevant Dirac phase in these works comes out to be different from that in our 
parametrization. 

Now we discuss unitarity triangles and neutrino oscillation formulae which get sim- 
plified using this parametrization. The unitarity triangles may be constructed using 
the orthogonality of different pairs of columns or rows of the mixing matrix. Infor- 
mation about the elements Ue2, f^es and f/^3 is obtained in the solar, reactor and 
atmospheric experiments and the most important unitarity triangles should include 
all these elements [23]. Two such unitarity triangles correspond to the orthogonality 
of second and third column z/2-'^3 and orthogonality of the first and second row l'e■'^^l■ 
The unitarity relation for the 1^2 m^ triangle is given by 

Ue2U:, + U,2U;, + Ur2U:, = (21) 

and the sides of this unitarity triangle to first order are given by 
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which satisfies Eq.(21) to first order. The invariant Jqp is 

JcP = Im{U,2U:,U;,U,s) = Im{Ur2U:,U:,U,s) = Im{U,2U;^U:^Urz). (23) 

The other unitarity triangle z/e.z/^ corresponds to the unitarity relation 

u,,u:, + u,2u:, + u,,u:, = o (24) 

and to first order the sides of this unitarity triangle are given by 
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which satisfy Eq.(24) to first order. The invariant Jqp is 

JcP = ImiU,sU:,U;,U,2) = ImiU,,U:,U;,U,s) = /m(f/^2f/:2f/;if/ei). (26) 

Now, we discuss the applications of this parametrization to neutrino oscillations. The 
probability of oscillation from fiavor z/^ to fiavor u/^, P{va -^ ^p) is given by 

i=l 

= 5,^ - 4 E Re{Ul,Up,U^,U;^) W A,, + 2 ^ Im{Ul,Up,U^,U;^) sin 2 A,, (27) 

where a, P = e, n, r, Ajj = (m^ — m'j)L/4E, L is the oscillation length and E is 
the beam energy of neutrinos. Expanding to second order in eij and A21 assuming 
A21 <^ 1 [23], we obtain the various vacuum oscillation probabilities. 
For e — )■ e, /i, r we obtain 
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The above equations give electron neutrino survival and disappearance probabilities 
to the second order in Cjj and A21. Note that these probabilities are independent of 
deviations from atmospheric mixing so that any deviation from maximal atmospheric 
mixing only appears at third order in these oscillation probabilities. 
For /i — )■ e, /i, r we obtain 
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New results have been announced by long baseline experiment T2K probing the 
i^fi -^ ^e appearance channel giving non zero reactor mixing angle [12]. The devi- 
ation from maximal atmospheric mixing does not appear upto second order in the 
oscillation probability -P(t'^ — > Ve)- Therefore, in this parametrization ^13 can have 
a large deviation from zero irrespective of the deviation of ^23 from ^ to a good ap- 
proximation. 
For r — )■ e, /i, r we have 
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P(^,^ ^ ^,^) = 1 _ sin^ A31 - ''^ '' + 46^3 cos^ sin^ Al,.{30) 

According to the above equations, some oscillation probabilities are identical upto 
the second order. This is the consequence of the simplifying assumptions A32 ~ A31. 
These oscillation probabilities will differ slightly when third order perturbations are 
considered. 



4 Summary 

Neutrino oscillation experiments suggest an atmospheric mixing angle very close to 
7r/4 and a small reactor mixing angle. We assume the lepton mixing matrix at zeroth 
order having ^23 = f and ^13 = 0. Deviations from maximal atmospheric mixing 
and vanishing reactor mixing are obtained through charged lepton corrections in 
terms of small perturbation parameters. Relatively large deviations for the reactor 
mixing angle from zero as indicated by T2K experiment can be obtained in this 
parametrization. In the zeroth order lepton mixing matrix, we keep the solar mixing 
angle general, so that the deviations from a particular mixing scheme e.g. TBM, 
GRl, GR2, HM and BM can be obtained by substituting the value of solar mixing 
angle. In this analysis, we have been able to reduce the number of complex phases 
to two by considering the Hermitian products of charged lepton and neutrino mass 
matrices, thus, resulting in considerable simplification of our analysis. The Jarlskog 
rephasing invariant measure of CP violation contains a single phase difference in the 
leading order which allows us to identify this phase difference with the Dirac-type 
CP violating phase in this parametrization. We have also obtained the formulae for 
the sides of unitarity triangles and vacuum oscillation probabilities. It is found that 
the deviation from maximal atmospheric mixing does not appear upto second order 
in the oscillation probability Pii^fi — )■ i^e) relevant for the measurement of reactor 
mixing angle. Therefore, in this parametrization ^13 can have a large deviation from 
zero irrespective of the deviation of 623 from maximality to a good approximation. 
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A Appendix 

Here, we list the second order corrections to the results given in the main text. The 
second order corrections to the first order mixing matrix elements are given by 
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The second order corrections to the mixing angles are given by 

A • n e23cos0(ei3 + ei2Cos(5i3 -0)) ei2sin^(5i3 -0)(l + 3e23COS( 
A sm fc'13 ^ 7= h 



Asin^ 



23 



Asin^ 



12 



72 ei32v^ 

^12 ^ ^13 + 4 ~ 2ei2ei3 cos(5i3 - 0) + £23 COS 2 

4^2 
3s'i2 - C12 cot 6';2 COS 2(^13 - 0)) - 2s[2eJ2 



e?3(csc^'i2 



+ 



\/2ci2ei2e23 COS0 - ei3(v/2ci2e23 cos(5i3 + s'^^eu cos((5i3 - 0)) 

2 
The second order correction to the Jarskog CP invariant is given by 
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The second order contributions to the z/2-^3 unitarity triangle are given by 
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